Generalized ghost dark energy in Brans-Dicke theory 
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It was argued that the vacuum energy of the Veneziano ghost field of QCD, in a time-dependent 
background, can be written in the general form, H-\-0{H^), where H is the Hubble parameter. Based 
on this, a phenomenological dark energy model whose energy density is of the form p = aH + j3II^ 
was recently proposed to explain the dark energy dominated universe. In this paper, we investigate 
this generalized ghost dark energy model in the setup of Brans-Dicke cosmology. We study the 
cosmological implications of this model. In particular, we obtain the equation of state and the 
deceleration parameters and a differential equation governing the evolution of this dark energy 
, model. It is shown that the equation of state parameter of the generalized ghost dark energy can 

cross the phantom line {wd = —1) in some range of the parameters spaces. 
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I. INTRODUCTION 



According to standard cosmology we believed that the universe is decelerating and one day it will stop expanding 
to reverse the process. The story has changed since 1998 with several exciting observations. The first significant 
result came from the distance measurement of the type la supernovaes [SNela], which states that the universe is 
^—1, experiencing a phase of acceleration |l|, 0| • This result has been confirmed repeatedly by several other observations 
^ ■ such as measurement of the anisotropies of the cosmic microwave background (CMB), spectrum by the Wilkinson 
y[ Microwave Anisotropy Probe (WMAP) Q, and by the measurement of the baryon acoustic oscillations (BAO) in the 
^ . Sloan Sky Digital Survey (SDSS) luminous galaxy sample 

In order to explain such an unwanted phase of acceleration in the Einstein's framework of gravity we need to 
assume that the universe is filled with an unknown type of energy which is called dark energy (DE) . This component 
CNJ ■ of energy have a negative equation of state parameter (EoS) w = ^ ~\ ^'^'^ derive such an acceleration. One 
main task for the theoretical physicists is to identify the nature of such DE. The simplest candidate is the famous 
cosmological constant with w = — 1, which is still one of the best models, agrees with observations. However, this 
candidate suffers the so called fine-tuning and the coincidence problems (5|. Further observations favor alternatives 
■ whose EoS parameter change with time. Simplest example of this class is the scalar fields which deserves time var yin g 
EoS parameters. An incomplete list of the scalar filed or time varying EoS parameter models can be found in |6l4l5l|. 
^ In addition to the DE approach, one can also explain the late time acceleration of the universe with modified gravity 

fS| ! [3 or inhomogeneous cosmology [l7j . In these approaches one assumes that the underlying theory of gravity should 
be modified in a way to lead a phase of acceleration at late time. 

Since any new model of DE has many unknown features and can lead new problems in the literature, our prior 
choice is to handle the DE problem without introducing new degrees of freedom beyond what are already known. 
Recently one class of such models has been attracted a lot of attentions, the so called "ghost dark energy" (GDE). In 
5-H ] this approach the Veneziano ghost field is responsible for the recent cosmic acceleration. The Veneziano ghost field 
was proposed to resolve the U(l) problem in QCD. The U(l) problem is that the Lagrangian of QCD has, in the 
massless limit, a global chiral U(l) symmetry, which does not seem to be reflected in the spectrum of light pseudoscalar 
mesons. The ghost field has no contribution in the vacuum energy in Minkowski spacetime but in the time dependent 
or non-flat backgrounds ghost filed has a non vanishing contribution to the vacuum energy proportional to Aq^piJ, 
where H is the Hubble parameter and Aq(^q is QCD mass scale In the GDE model the vacuum energy of the 
ghost field can be taken as a dynamical cosmological constant [3, [13 ■ Different features of GDE have been explored 
in ample details [20l - |23j . In [23], the author discussed that the contribution of the Veneziano QCD ghost field to 
the vacuum energy is not exactly of order H and a subleading term appears due to the fact that the vacuum 



expectation value of the energy- momentum tensor is conserved in isolation 25| . They argued that the vacuum energy 
of the ghost field can be written as H + 0{H^), where the subleading term in the GDE model might play a crucial 
role in the early evolution of the universe, acting as the early DE [26l|. This term can also lead to a better agreement 
with observations ,27^ . 
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On the other hand, in recent years scalar tensor theories have been reconsidered extensively. The reason comes from 
the fact that scalar fields appear in different branches of theoretical physics as a consistency condition. For example, 
in the inflation theory proposed by Alen Guth a scalar field was introduced as an inflator [i^, or the low energy 
limit of the string theory leads to introducing a scalar degree of freedom. One important example of the scalar tensor 
theories is the Brans-Dicke (BD) theory of gravity which was presented by Brans and Dicke in 1961 to incorporate 
the Mach's principle in the Einstein's theory of gravity [1^. This theory also passed the observational tests in the 
solar system domainjs^]. It was shown that some features of original GDE in BD cosmology differ from Einstein's 
gra vity. For example while the original GDE is instable in all range of the parameter spaces in standard cosmology 
[23} , it leads to a stable phase in BD theory ^l|. All the above reasons motivate us to investigate the generalized 
GDE (GGDE) model with subleading term in the framework of BD theory. In the next section, we review the 
GGDE model in standard cosmology. In section Hill we generalize the study to the BD framework. We finish with 
closing remarks in section IIVI 



II. GENERALIZED GHOST DARK ENERGY MODEL 



For flat FRW universe filled with GDE and pressureless dark matter, the first Friedmann equation may be written 

as 

H^ = ^{pm + Pd) (1) 

where pm and po are, respectively, the energy densities of pressureless matter and ghost dark energy. The generalized 
ghost energy density is [26i] 

PD=aH + I3H^ (2) 

where a is a constant of with dimension [energy]^ , roughly of order Aq(^q and Aqcd ~ \QQMeV is QCD mass scale 
and /? is another constant with dimension [energy]^ . The motivations to consider this form of energy density is that 
this form of varying cosmological constant could be a possible candidate to solve the two fundamental cosmological 
puzzles, the fine tuning and the coincidence problems [32| . 
We define the dimensionless density parameters as usual, 

o n PD a + fiH 

Per Per 6M^H 

where the critical energy density is pcr = 'SH^Mp. Thus, the Friedmann equation can be rewritten as ilm + — 1- 
The conservation equations read 

p,n + iHp„, = 0, (4) 
PD + iHpD{l + WD) = 0. (5) 

Integrating Eq. (|H) , we find 

Pm = /Omo(l + 2:)^, (6) 

where z = l/a— lis the redshift function. Thus, Friedmann equation ([T]) can be written 

H^T - aH ^ p,r,o{l + zf (7) 
where F — 3Mp — /3. Solving the above equation we find 

H{z) = ^ (a ± Va'+4rp™o(l + ^)3) (8) 
Using the fact that pmo — SMpHQ^lmQ = (F + 13)11^0,^0, the above equation can be further rewritten as 

H{z) = ^ (a ± ^a^+^T{T + l3)H^a^o{^ + zf^ (9) 
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Using the energy density ratio relation u = Pm/ Pd, the Friedmann equation can be also written as 

SM^iJ^ = {l + u)pD. 



(10) 



Taking the time derivative of relation ([2]) and using the Friedmann equation we find pu = (a + 2j3H)H. Also, from 
the Friedmann equation as well as continuity equations we have 



H 



2M|> 



{1 + U + wd)p_ 



D 



Substituting po into Eq. ([5]), after some simplifications, we find the EoS parameter of the GGDE as 

1 

WD 



where ^ = — . It is easy to see that at the early time where VtD we have 



WD 



(11) 



(12) 



(13) 



while at the late time where rjjj — )• 1 the generalized ghost dark energy mimics a cosmological constant, namely 
wd = — 1- This behavior is similar to the original GDE [21|. When ~ 13 ~ 0, one recovers the EoS parameter of the 
original GDE ^211] 



WD 



2-n, 



(14) 



Also, the deceleration parameter is obtained as 



H 



-1 + -nD{i + u + ujD) 



(15) 



where we have used Eq. ([TT|) and relation pD = 3MpH^. Substituting wd from p^ . we can further simplify q as 



q^^-^nD{2-nD + 2^H (1 - ^d))-^ 



(16) 



At the late time where the dark energy dominates {^d 1) we have q — —\. Taking the time derivative of in 
Eq. 0, we find 



n, 



f a \ H aflDH 



Using relation flD ~ ^Jina^ obtain 

d^l D ^D 



din a 



S f Qd 



\3Mj,J a + l3H 



(1 + 9) 



(1 + 9) 



2\1+^H 



D 



{2-nD) + 2^H{i~nD) 



(17) 

(18) 
(19) 



This is the equation of motion governing the evolution of GGDE. The evolution of H in Eqs. and ([T5| is given 
byEq. (HI). 



III. GENERALIZED GHOST DARK ENERGY IN BD COSMOLOGY 



The action of BD theory, in the canonical form, can be written [33 



(20) 
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where R is the scalar curvature and is the BD scalar field. Varying the above action with respect to FRW metric, 
f/^ii, and the BD scalar field (j), we obtain the following field equations 

+ + + (21) 

4^*' (4 + 4) - - i** - K' + ^) ^ 

i + 3H4,-^(- + H' + -^),P = 0. (23) 

In the following, we consider the flat FRW universe, thus we set k = 0. At this point the system of our equations is 
not closed and we still have a freedom to choice the scalar field. In the framework of BD cosmology the BD scalar 
field (j) is usually assumed to has a power law relation in terms of scale factor, namely, 

= Mty. (24) 

A case of particular interest is that when e is small whereas uj is high so that the product euj results of order unity 



3J|. This is interesting because local astronomical experiments set a very high lower bound on uj [35]; in particular, 
the Cassini experiment implies that uj > 10^ ^30., 36]. Taking the time derivative of relation ([M)) . we obtain 



e- = eH. (25) 
a 



Using Eqs. ([M)) and (pS)) . the first Friedmann equation (|21l) can be rewritten [k = 0) 



H'{l-—e' + 2e) = —{pD + Pm). (26) 



3 '30- 

We introduce the fractional energy densities corresponding to each energy component as 



where we have defined the critical energy density as 

The reason for this definition comes from the fact that in BD theory, the non-minimal coupling term (jP'R replaces with 
the Einstein-Hilbert term R/G in such a way that — 27r0^/aj, where Gog is the effective gravitational constant 
as long as the dynamical scalar field varies slowly. 

Using generalized ghost energy density ^ we can rewrite Eq. (^5)) as 



^D^:^ia + I3H), (30) 



4a; 

Using definitions (P7)) and (pS)). Eq. ([^ can be expressed as 

no + ^m^i (31) 

where we have defined 

7=l-|^e2+2e (32) 

Clearly for £ = (a; — oo) we have 7 = 1. In this case the BD scalar field becomes constant and Einstein gravity is 
restored. Taking the time derivative of Friedmann equation (j26p we find 



2HH-f=^ipD+Pm)-^{pD+Pm). (33) 
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Using the continuity equations, as well as Eq. (1251) we get 



2Hi^-^{1 + u + lod)-^^Pd{1 + u). (34) 



Dividing by H"^, we have 



Using the fact that (1 + u)VId — -f we obtain 



-SnD{l + u + ujD) -^enoil + u) . (35) 



H / 3 3 ^dujd\ 



Finally for the time derivative of ghost energy density, we obtain 



PD = ia + 2I3H) H^-H^e + ^ + (a + 2/3iJ) . (37) 



Inserting this relation in Eq. ([5]), it is a matter of calculation to show that 

WD = — — . (38) 

(27-f7D) + 2M(7-Op) 

When 13 = one recovers 



= TTT-^TTT - 1 ) • (39) 



{2^ -no) V3 

which is the EoS parameter of the original GDE in BD theory presented in [22]. On the other hand, in the absence 
of BD scalar field e = (7 = 1) we obtain the result of the previous section, namely 

^ i2-nn)+'^ii-n,y ^''^ 

The solar-system experiments give the lower bound for the value of w to be a; > 40000 [30*1. However, when probing 
the larger scales, the limit obtained will be weaker than this result. It was shown ^] that w can be smaller than 
40000 on the cosmological scales. Also, Wu and Chen [s^] obtained the observational constraints on BD model in a 
fiat universe with cosmological constant and cold dark matter using the latest WMAP and SDSS data. They found 
that within 2a range, the value of ui satisfies ui < —120.0 or w > 97.8 [s^]- They also obtained the constraint on the 
rate of change of G at present 

G 

- 1.75 X IQ-^^yr-^ < n < 1-05 x 10"^2yr"\ (41) 



at 2(7 confidence level. As a result in our case with assumption ([24| we get 

^ = - =ei/ < 1.05 X 10"i2y^-i^ (-42) 
G (j) 

This relation can be used to put an upper bound on e. Assuming the present value of the Hubble parameter to be 
Hq ~ 0.7, we obtain 

e < 0.01. (43) 

The GGDE model in BD framework has an interesting feature compared to the GDE model in Einstein's gravity. It 
was shown that in standard cosmology based on Einstein's theory, the EoS parameter of the noninteracting GGDE 
cannot cross the phantom line wd = —1 and at the late time where ^Id 1 approaches — l[2lj. Choosing fi^io = 0.72, 
Hq = 0.7 for the present time and a suitable choice of ^, this inequality valid provided we take e < 0.01 in a narrow 
range which is consistent with observations. This indicates that one can generate a phantom-like EoS for the GGDE 
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in the BD framework. One should note that increasing ^ can exclude crossing the phantom line which indicates a 
negative contribution of the subleading term with respect to the leading term H in the energy density. 

Since the dynamic of the universe should be discussed in term of the effective EoS parameter thus in addition to 
the EoS parameter of the GGDE we also study the effective EoS parameter, WcS, which is defined as 

w., = ^ = (44) 

Pt PD + PM 

where pt and Pt are, respectively, the total energy density and total pressure of the universe. As usual we assumed 
the dark matter is in the form of pressureless fluid (P/\/=0). Using relation for the flat case one can find 



^<=ff = '^D = 7f- Lh , r. . (45) 



7 (27 - f^i,) + ^ (7 - f^D) 

It is also interesting to study the behavior of the deceleration parameter defined as 



Substituting Eq. ([TT|) in the above relation one can easily reach 

1 3 flnWD 



2 27 



(47) 



Inserting Eq. ^ into (gT]) yields 



When /3 = we obtain ^ 

In the limiting case e = (w — >■ 00) we have 7=1 and hence the BD scalar field becomes trivial; as a result Eq. pS)) 
reduces to its respective expression in flat standard cosmology obtained in the previous section 

q= 1-^^0(2 -no +'^{1-^0)) (50) 
2 2 \ a J 

Let us study some special cases of interest for the deceleration parameter q. If we take ri^o = 0.72, Hq ~ 0.7 and for 
the present time and choosing e — 0.002,^ = 0.1 and w = 10000 we obtain go = —0.35, which is consistent with the 
present value of the deceleration parameter obtained in (38j . Transition from deceleration to acceleration occurs at 
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fiu = 0.46. It is worth noting that GGDE results in a smaher rate of acceleration in comparison with the GDE then 
this models leads a delay in different epoches of the cosmic evolution with respect to the original GDE (for instance 
the GDE with choice of a same set of parameters lead qo = —0.37). For a better insight about these two models we 
plotted Wef f and q for both of models in Fi^ 

Finally, we obtain a differential equation governing the evolution of GGDE from early deceleration to the late time 
acceleration. Following the method of the previous section we find 

^^^^{q + l)-2nDeH (51) 

= ^dh( ^(g+l)-2e). (52) 



Inserting q from (HS)) and using relation (Id = Hil'j^, we obtain 

^'d = (^^^ (1 + 9) - , (53) 

where the prime denotes derivative with respect to x = In a and q is given by Eq. (I48p . In the limiting case (3 = 
one recovers the result obtained in (2^] . 



IV. CLOSING REMARKS 



A phenomenological GGDE model whose energy density is of the form p = aH + f3H^ was recently proposed to 
explain the observed acceleration of the universe expansion. This model originates from the fact that the vacuum 
energy of the Veneziano ghost field in QCD is of the form, H + 0{H^). It was shown that the difference between 
vacuum energy of quantum fields in Minkowski space and in FRW universe can play the role of observed dark energy. 

In this paper we first studied the cosmological implications of the GGDE model in standard cosmology. We found 
that in this model, the universe approaches to a de Sitter phase at late times. Then, we extended our study to the 
BD cosmology. We can classify our achievements in two categories. The first is that in the BD framework the GGDE 
can cross the phantom line with suitable choice of free parameters while in the standard cosmology we found a de 
Sitter phase as the fate of the universe. The second result is that the subleading term can lead a delay in different 
epoches of the cosmic evolution. We also discussed this result explicitly by a numerical evaluation and showed that 
taking a same set of parameters result qo = —0.35 for the GGDE while for GDE gives qo = —0.37. It is easily seen 
from Fig IT] that GDE (dashed curve) enters the acceleration phase sooner than the GGDE (solid curve). This point 
was also addressed in standard cosmology f26l| as the negative contribution of the subleading term in the energy 
density. 
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